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Abstract 

We compute the dispersion laws of chaotic periodic systems using the semiclassical 
periodic orbit theory to approximate the trace of the powers of the evolution operator. 
Aside from the usual real trajectories, we also include complex orbits. These turn out 
to be fundamental for a proper description of the band structure since they incorporate 
conduction processes through tunneling mechanisms. The results obtained, illustrated 
with the kicked-Harper model, are in excellent agreement with numerical simulations, 
even in the extreme quantum regime. 
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As is well known, the quantum motion of electrons in periodic potentials gives rise to en- 
ergy bands. A precise determination of the density of states is necessary in order to obtain 
an appropriate description of the elementary physical properties, like the transport and mag- 
netic properties. However, a direct calculation of the dispersion laws is in general difficult, in 
particular if the underlying classical dynamics is chaotic. 

The purpose of this letter is to consider the band structure of classically chaotic periodic 
systems from the point of view of periodic orbit theory. In this approach, the trace of the evolu- 
tion operator - which can be related to the density of states - is expressed, in the semiclassical 
limit, as a sum over classical periodic paths 0,0- As shown below, the dependence of the dis- 
persion laws on the quasimomenta (or Bloch angles) is related to the winding numbers of the 
periodic orbits (p.o.'s) around the elementary cell. Trajectories with winding numbers different 
from zero are associated to open diffusive processes through the periodic lattice. However, in 
some cases these processes are classically forbidden, or not well described by the real p.o.'s. 
The relevant features of the band structure can however be recovered going one step beyond 
this approximation by including some tunneling effects in the semiclassical trace formula. This 
is done by taking into account the complex p.o.'s (by complex or ghost p.o.'s we mean the pe- 
riodic solutions of Hamilton's equations having at least one complex phase-space coordinate). 
Moreover, the inclusion of ghost orbits extends the range of validity of the semiclassical ap- 
proximation. As an illustration, at the end of the Letter we will present an example in which, 
in spite of being far from a classical behaviour, we are able to compute the (quasi-)energies as 
a function of the Bloch-numbers (and hence the transport properties of the system) in terms of 
complex p.o.'s, while only taking into account the real p.o.'s gives a very crude approximation 
to the true solution. 

In the standard WKB theory - valid for integrable systems - complex classical paths have 
been used in the past in different contexts, as in the study of the band structure of one- 
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dimensional periodic potentials. To our knowledge, the existence and relevance of ghost orbits 
in the semiclassical trace formula was pointed out for the first time in |2j. More recently, 
these orbits were studied in more detail in Refs. 0, [5] in connection to the trace of the evolution 
operator of the kicked top and the standard map, respectively. Here, we also consider a dynam- 
ics described by a discrete (chaotic) map, inspired by the study of the motion of an electron 
in a two-dimensional periodic potential in the presence of a uniform and constant magnetic 
field perpendicular to the crystal plane. For a strong magnetic field, the lowest Landau level 
approximation leads to the motion of a ID-particle whose phase space is a 2D-torus T. If 
Planck's constant is rational, h = M/N, then the corresponding evolution operator U is doubly 
periodic and its spectrum is made of N bands e a ( 6) given by U\ip a ( 6) >= exp(ie a ( 9))\ijj a ( 6) >, 
a = 1, . . . , N. For simplicity, in the following we take M = 1 and T will be of unit length 
in both directions. The Bloch parameters 9 = (61,82) are associated to generalized periodic 
boundary conditions ij)(q + 1) = exp (\9\) i/j(q), ij){p — 1) = exp (i# 2 ) i>{p)- The classical limit is 
obtained when N 00. 

In order to mimic the motion of the electron as the magnetic field is lowered, we consider 
the kicked Hamiltonian H = f(p) + g(q) $(t ~ w), which leads to a discrete map on the 

n=— 00 

torus 

ft+i = Qi + rf'iPi+i) modi 
Pi+i = Pi ~ rg'(qi) modi, 
where (qj,Pj) G [0, 1[ 2 are the phase-space coordinates of the particle at time t = jr while / 
and g are two smooth functions of period one. One can interpret the two integers 

= h + r/W)] > =~\Pi- r 9 '(q i )} 

(square brackets denote integer part) as the winding numbers of the orbit around the torus for 
the iteration (qi,Pi) — > (qi + xiPi+i) m the q and p directions, respectively. As we shall see in 
the following these two integers - which arise from the toroidal topology of phase space by the 
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identification of points differing by an integral number of phase-space cells - play a crucial role 
in the semiclassical interpretation of the bands (see also Ref.||). Quantum mechanically the 
map is implemented by the one-step evolution operator 

JJ = e -iTf(p)/h e -iTg(q)/h^ ^) 

In order to make explicit the dependence on the parameters 9, instead of U we consider 
the operator U(9) which is obtained from Eq.(§) by making the substitutions q — > q + h9z, 
p —> p + H9i 0. Then the dispersion laws e a (9) are given by the roots z a (9) = exp(ie a (9)) 
a — 1 , . . . , N of the spectral determinant 

N 

det(z-U{9))=J2a k {9)z k , (3) 

k=0 

where a^(9) = 1 by construction. The coefficients of the characteristic polynomial (^) can be 
obtained from the trace of the first N powers of U through the formula 

1 k 

a N _ k = -- Y,CLN-k+ntr{U n ), k = l,...,N. (4) 

K n=l 

The interest of (|]) comes from the fact that in the semiclassical limit iV — > oo the trace 
of the powers of U can be computed in terms of a finite number of classical p.o.'s (which are 
assumed to be isolated) [|I|, [|. In our case, by considering the special topology of phase space, 
we get 

tr (U n ) =i n J2 T =L exp (~S + iw p 9 2 + iw q 9 x - i|z/) . (5) 



p.o. y/\ det(M - 1) 

In this formula the sum extends over all the real p.o.'s of the map ([IJ) whose period T is an integer 

n-l n-1 

fraction of n, w p and w q are the total winding numbers of the p.o., w p = J2 W J , w q — J2 , 

3=0 ' 3=0 

S is the classical action of the p.o. 

n-l 

S = {~ T (9(Qj) + fiPj)) + Pj+i{Qj+i ~ Qj) ~ Wp )( lj - w^Vj) , 



3=0 



j=n ^ -Tg"(qj-i) 1 

is the monodromy matrix (related to its stability) and v is the number of negative eigenvalues 
of the 2n x 2n matrix (assumed to be non singular): 

/ ( d 2 s \ ( d 2 s \ \ 

\dPidPj J i,j=l,..., n \ d Pi d H h,j=l,...,n 

( a 2 s \ ( d 2 s \ 

\ \dqidpj J it j=i t ... t n \9^9<lj Ji,j=l,...,n J 

The information concerning the band structure (i.e., the dependence on 9 of the dispersion 
laws) is thus recovered, semiclassically, through the factor exp (i {w q 9\ + Wp^))- This factor, 
which is different from zero for open diffusive p.o.'s, can be interpreted for each 9 as a unitary 
representation of the homotopy group of T. The fact that the representation of that group 
appears in tr (U n ) when one deals with path integrals (or their discrete equivalent) in multiply- 
connected spaces is well known 0. 

Because U is a unitary operator, the coefficients in Eq.(Rf) satisfy the self-inversive property 
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a N - k = exp (i6) a k , k = 0, . . . , N (6) 

where 6 = ttN + J2^a(9). This property, together with Eqs.@-@, allows to compute the 
band spectrum in terms of the classical p.o.'s of length less than [N/2]. 

Eq.(|5]) was obtained from the exact expression of tr (U n ) replacing the sums by integrals 
using the Poisson summation formula and then computing the integrals by the stationary phase 
approximation. One step further is to deform the contour of integration of the path integral in 
order to reach the stationary orbits lying in the complex plane. The outcome of this procedure 
is to include in the summation flj|) complex p.o.'s, whose contribution is exponentially small. 
However, as we shall now see, they cannot be neglected if H is not too small or if we are close 
to a (first order) bifurcation of a p.o. . 
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For the sake of definiteness, we consider the kicked-Harper model 

g(x) = f(x) = -7cos(27rx)/(27rr). (7) 

Phase-space plots of the classical trajectories of the kicked-Harper map for different values of 
7 > can be found in |7j]: for 7 — > + the dynamics is integrable (it tends to the Harper 
Hamiltonian); for 7 ~ 0.4 it is a mixture of regular and chaotic orbits, while for 7 > 0.6 
unstable chaotic trajectories dominate. At the quantum level, some spectral and transport 
properties of the model were studied in |f7|, jTT | . 



We now compute the trace of U for this model. The exact result is 

tr U = N { J o ( 7 A0 + 2 f] e- iwN ^ 2 C os(w6 1 )J wN ( 1 N)) x {9 l # 6 2 } , (8) 



W = l 



where the second curly brackets in the r.h.s. is identical to the first one except for the fact that 
61 is replaced by 62- In order to compute (§) semiclassically, we need to find all the p.o.'s of 
period one of the map. Their location (q,p) G [0, 1[ 2 is given, from (|l|) and (0), by the solutions 
of the set of equations 

sin(2vrg) = w P h 
sin(27rp) = Wg/j, 

where (w q ,w p ) G 2L 2 . Let us briefly discuss the different solutions of Eq.flj^). Assume for the 
moment that both w q and w p are positive. A set of winding numbers (w q ,w p ) corresponds 
to a real p.o. if w q ^ p < 7. In that case there are four different solutions of Eq.@ in [0, 1[ 2 : 
{Qo, Po), (go, V 2 -Po), (I/ 2 - %iPo) and (1/2 - q , 1/2 -p ), where q = arcsin(wp/7)/2vr and 
p = arcsin(t(7 (? /7)/27r. On the other hand, (w q ,w p ) correspond to a complex p.o. if at least 
one of its winding numbers is larger than 7. If both integers are larger than 7, there are four 
different p.o.'s located at 

a = l/4±iArchr?y;„/V)/27r 

(10) 



p = 1/4 ± iArch (wp/j) /2vr; 
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while if Wi < 7 < Wj one of the coordinates is real and the other one is complex. 

For negative values of w q and/or w p , the same kind of analysis can be done. The global 
description of the period-one p.o.'s is therefore as follows. For a fixed 7 > 0, there is an 
infinite number of complex p.o.'s (all integers (w q ,w p ) whose modulus is larger than 7). Since 
the imaginary part of the coordinate is proportional to Arch |i«i/7|, the complex p.o.'s with 
large winding number \wi\ » 7 are deeply located in the complex plane. As 7 increases 
(starting from 1/2, for example), a complex orbit becomes real if 7 becomes larger than the 
modulus of both winding numbers. This will happen at each integer value of 7 through a first- 
order bifurcation (fold catastrophe), where several orbits coalesce in at least one phase-space 
coordinate. The total number of real p.o.'s increases with 7 as (4 [7] + 2) 2 . 

Concerning the stability of those orbits, generically trM = 2 at 7„ = n G IN* for all the 
real p.o.'s emerging at that value of 7. When 7 is increased, half of them are initially stable, 



the other half unstable; but very rapidly they become all unstable JT2[ . 

Using the information described above we have computed semi-classically the trace of U 
(for fixed 7 and arbitrary N) including all the real as well as complex p.o.'s. The result is 

trU = —F(9 1 )F(9 2 ) + 0(1/N) (11) 

7T7 



where 

F{9) 



w=0 \/COS[X v 

00 



clw Nn/2 COs(w6>) e -Ar( Mfe - 7S inh (y w )) ^g) 

w=[H+i Jsinh (y w ) 



(13) 



and 

x w = arcsin(w/7) 
y w = Arch(w/7). 

In Eq.flnp, F(9 2 ) is related to the g-coordinate of the orbit, while F{9{) contains the information 
about its momentum. The first sum in Eq.(^) (w < [7]) comes from p.o.'s having a real 



coordinate (which reproduces the oscillating part of Bessel functions occurring in Eq.(|3|)), while 
the second sum (w > [7]) comes from those having a complex coordinate (reproducing the 
exponentially decreasing part of the Bessel functions) . The symbol indicates that the term 
w = must be multiplied by a factor 1/2. The product of both functions in Eq . (|TT|) gives 
rise to four different sums, one related to real p.o.'s, the other three to complex ones. The 
dependence on 9 in Eq.(|TI|) is the same for both kinds of orbits, but the complex orbits are 
exponentially dumped by a factor (or a product of factors) of the form exp(— NSj), where 
Sj = wy w — 7sinh (y w ) > 0. Because of that, complex p.o.'s with large winding number can in 
general be neglected. On the other hand, as 7 approaches some integer value n from below, the 
contribution of the complex p.o.'s becoming real at 7„ = n > is particularly important, since 
the imaginary part of their action tends to zero. In fact, they remain important for parameter 
values relatively far from j n , since Sj ~ \j2/n(j n — 7) 3//2 /37r as 7 — > 7„ from below (the exponent 
3/2 was also found for the kicked top and the standard map [§], Q). Moreover, because of the 
1 / \fn dependence, the parameter interval in which complex orbits are important increases with 
7 (at a fix N). Exactly at 7 = 7„ Eq. flTTD is not defined since the denominator of both sums in 
F{9) vanishes due to the coalescence of several orbits. In order to avoid divergences we must 
improve the approximation, and include third order terms in the computation of the integrals 
for 7 ~ 7„ using Airy functions [13|, U . Including these corrections, we find that (|Tl|) reproduces 
extremely well the exact tr U, even for small values of N. 

In order to illustrate this, let us consider the extreme quantum limit N = 2. For the kicked- 
Harper map, it can be shown that det U = 1 for arbitrary 7, N and 9, and then = tcN 
in Eq.(^). From this it follows that for N = 2 all the information concerning the spectrum 
of U is contained in tr U, since Eq. ([|) reduces to a second-degree polynomial with coefficient 
ct2 = 1, Oi = — trlf, (Zq = 1. The two dispersion relations obtained from that equation are 
e ±($) — =t arccos(tr U( 9)/2). Although our results are valid for arbitrary 7 > 0, the role of 
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the complex p.o.'s is particularly stressed in < 7 < 1. In that parameter interval there 
exist only four real p.o.'s of period one, located at (0,0), (0,1/2), (1/2,0), (1/2,1/2). But 
these orbits do not introduce any 9 dependence in txU (and therefore in e±(8)) since for all 
of them w q = w p = 0, leading to flat bands. In fact, the exact dispersion laws for N = 2 have 
a strong and non-trivial ^-dependence; in particular, for an arbitrary 7 there is at least one 
linear contact (diabolical point) between the two bands. In order to introduce semiclassically a 
6 dependence in the dispersion laws one must include complex p.o.'s. Fig. 1(a) plots the exact 
result and the semiclassical dispersion laws computed from Eq.([0]) with and without complex 
orbits for 7 = 0.86 and 82/271 = 0.5 as a function of 9i/2ti (only complex p.o.'s up to winding 
number two were included). The result obtained with complex p.o.'s is extremely good even 
near the diabolical point, which is the worst situation for semiclassical analysis. Part (b) of 
that Figure shows the case 7 = 4.86. 

It is also possible to compute, using only the trf/, the dispersion laws for N = 3 since 
by the self-inversive symmetry and the fact that O = 3n the coefficients of the characteristic 



polynomial are now given by 03 = 1, 02 = —trU, a\ = trU, a = 1. Fig.l(c)-(f) illustrate the 
results obtained in this case for the dispersion laws and the trace of U. 

As a final remark, let us point out that the approach presented in this Letter is also relevant 
in systems where the topology producing the band structure is not toroidal, but for example 
cylindrical. This occurs in particular in the context of persistent currents, where the Aharonov- 
Bohm flux threading the ring plays the role of the Bloch-parameter and the winding number 
of the p.o.'s is precisely the winding number around that flux (see Ref.[[H|] for a more detailed 
discussion) . 
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FIGURE CAPTION 



FIG. 1: Band structure and trace of the evolution operator for the kicked-Harper model, (a) 
dispersion laws e( 6) as a function of 9i/2n for 62/271 = 0.5, N = 2 and 7 = 0.86. Full 
line: exact; dotted line: semiclassical theory including the real and complex p.o.'s; dashed 
line: semiclassical theory including only the real p.o.'s. (b) same as in part (a) but for 
7 = 4.86 and 9 2 /2n = 0.25. (c)-(d) same as in part (a) but for N = 3, 6 2 /2tt = 0.25, 
7 = 2.86 and 7 = 52.1, respectively, (e)-(f) real and imaginary part of tr (U) as a function 
of 7 for 6i/2n = 0.55 , 9 2 /2n = 0.32 and N = 3. The semiclassical result using only the 
real p.o.'s is not shown in (e)-(f). Except near the diabolical point in part (a), in all the 
other plots the dotted line cannot in general be distinguished from the full line. 
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